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Traditionally, since the discovery of the Josephson effect in 1962, the magnetic diffraction pattern 
of planar Josephson tunnel junctions has been recorded with the field applied in the plane of the 
junction. Here we discuss the static junction properties in a transverse magnetic field where de- 
magnetization effects imposed by the tunnel barrier and electrodes geometry are important. Mea- 
surements of the junction critical current versus magnetic field in planar Nb-based high-quality 
junctions with different geometry, size and critical current density show that it is advantageous to 
use a transverse magnetic field rather than an in-plane field. The conditions under which this occurs 
are discussed. 
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INTRODUCTION 

It is well known that a magnetic field H modulates the 
critical current Ic of a Josephson Tunnel Junction (JTJ) 
[l| . Indeed the occurrence of such diffraction phenomena 
Ic{H) is one of the most striking behaviors of JTJs|2]. In 
the case of a planar JTJ, it was Josephson himself [3] who 
pointed out that the gradient of the Josephson phase (j), 
which is the difference between the complex wavefunction 
phases in the electrodes, can be expressed as 



insensitive to transverse fields, this is the most obvious 
choice of the magnetic field orientation. 

For the reasons above, a magnetic field parallel to the 
barrier of planar JTJs is applied in the practical applica- 
tions of the Josephson magnetic diffraction phenomena 
such as, for example, the suppression of the d.c. Joseph- 
son effect in SIS mixers for photon detectionj^ and in 
specially shaped JTJs for particle detection^], the mag- 
netic biasing of a flux flow oscillator the tuning of 
resonant fluxon oscillators [ij. 
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where n is a unit vector normal to the insulating barrier 
separating the two superconducting electrodes, fiQ is the 
vacuum permeability and $o = h/2e is the magnetic flux 
quantum. If the two superconducting films have thick- 
nesses o?i,2 and London penetration depths Xli,2 and tj 
is the barrier thickness, then the effective magnetic pen- 
etration de is given byj; 



de = tj + Ail tanh h tanh ■ 
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which, in the case of thick superconducting films {di >> 
Xli), reduces to de « Xli + Xl2- 

bmce RowellfHi in 1963 made the first experimental veri- 
fication of EqlJl a large number of theoretical and ex- 
perimental papers have been devoted to the study of 
magnetic diffraction patterns, in various Josephson junc- 
tions. Nowadays the Ic{H) curves for planar JTJs hav- 
ing the most common geometrical and electrical param- 
eters is fully understood (see, for example, Cpt.4 and 5 
of Ref. [H). It is important to point out that nearly all 
work was done with the external magnetic field applied 
in the barrier plane. In fact, since Eq[T] states that (p is 



TRANSVERSE MAGNETIC FIELD 

In 1975 Rosenstein and Chen[lo| first reported on the 
effect of a transverse magnetic field on the critical current 
of a JTJ with overlap geometry. Among other things, 
they showed that the value of the junction critical field 
He at which the magnetic diffraction pattern first goes 
to zero, changes with the inclination of the field with re- 
spect to the barrier plane, the minimum being obtained 
when the field is transverse. This was the first experi- 
mental observation that transverse fields could be more 
efficient that parallel ones in modulating the Josephson 
current. Soon after, Hebard and Fulton'll"| correctly in- 
terpreted the findings of Ref. Jil] in terms of stationary 
screening currents which develop when a superconductor 
is subjected to an external magnetic field. To better un- 
derstand the mechanism through which also a transverse 
field is able to modulate the critical current of a planar 
JTJ, let us consider first a single isolated superconduct- 
ing film immersed in a uniform static magnetic field H± 
perpendicular to its surface. This system has received a 
continuous interest over the years and here we only recall 
the main features. For a deep treatment of this topic we 
remand to Ref. I4I and references therein. We assume 
that the film thickness d is larger than its London pene- 
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tration depth Al and that the field everywhere is much 
smaller than the critical field which would force the film 
into the intermediate or normal state, i.e., that the film is 
in the flux-fi-ee Meissner regime. At the top and bottom 
film surfaces, the flux lines are excluded from the interior 
of the film where H = 0. In fact, due to the screening 
currents Js(= VxH), they bend as they approach the 
film surface, flow along the film surfaces, concentrate at 
the film edges, and bent backward. Due to continuity, 
Hn, the component of H normal to the surface, may be 
taken to be zero, while its tangential component Ht de- 
cays exponentially inside the film on the scale of A^. 
The knowledge of the distribution of the magnetic field 
lines around the film requires a self-consistent solution of 
a magnetostatic problem combining the London equation 
(H + AlV^H = 0) in the superconducting film and the 
fourth Maxwell equation (V x H = 0) in the empty space 
around the film with boundary conditions appropriate to 
the film surface geometry. This problem can be solved 
analytically only for simple axially-symmetric cases such 
as, for example, that of an ellipsoid of revolution with the 
axis of revolution parallel to the applied field H±^ [l^ . 
If the film width w is much less than its length, but 
much greater than its thickness d, then we can approx- 
imate the film as a elliptical oblate cylinder of infinite 
length whose cross-section has axes w and d\ with the 
applied magnetic field H± directed along the minor axis, 
then i?tj the component of H tangent to the surface, 
only depends on the angle (3 with respect to the mi- 
nor axis: Ht/H± = 1/^cos^ /3 -I- (d/w)^ sin^ /3, whose 
maximum value w/d » 1 occurs at the cylinder edges 
(P = ±7r/2)[l3|. If the film width w is comparable to 
its length, the film can be approximated by a disk whose 
diameter w is much greater than its thickness d, with the 
external field applied parallel to its axis. In this case, for 
symmetry reasons, on the disk top and bottom surfaces 
the magnetic field lines are radial and Ht only depends 
on the distance r from the disk center; it is null at the 
center of the disk and increases as we move outwardfisj. 
The surface or sheet current density js, defined as the 
screening current density Jg integrated over the speci- 
men thickness, equals in magnitude Ht and is everywhere 
orthogonal to H. Numerical simulations carried out for 
a Nb disk having Xl = 90nm show that the shape of 
Ht{r) only depends on the disk aspect ratio w/d, as far 
as d > > Ai, . We found that at the disk border Ht is sev- 
eral times larger than the applied field i?_L, as shown in 
Fig.l for three values of the disk aspect ratio w/d = 10, 
100 and 1000, with d = 1/im. Both film approximations 
lead to conclude that a thin superconducting film of any 
geometry in a transverse field produce a magnetic field: 
i) whose orientation on the film end surfaces is paral- 
lel to surfaces themselves; ii) whose direction near the 
borders of the film end surfaces is perpendicular to the 
borders themselves; iii) whose strength is proportional 
to the transverse applied field intensity and exceeds its 
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FIG. 1: Computed normalized magnetic induction field 
Ht / at the border of a Nb superconducting disk in a per- 
pendicular field //x. The disk thickness is d = 1/im and 
diameters w = 10, 100 and IQOO^tm. 2r/w is the normalized 
disk radius. 



value near the borders of a film with a large aspect ratio 
w/d. 

Now, we can consider the situation in which two super- 
conducting films partially overlap to form a planar JTJ. 
If the tunnel barrier is very transparent, then the screen- 
ing currents cross the barrier and the two films act as 
one single fused fllm. In the opposite case the screening 
currents in the two fllms are independent on each other 
and no cross talk is allowed. In the intermediate situ- 
ations, a fraction of the screening currents in one film 
crosses the barrier and circulates in the other film and 
viceversa. It is clear that, for a given transverse field and 
a given junction and electrodes geometry, the transition 
from the fused to the independent fllms regime can be 
controlled via the barrier transparency, i.e., in our case 
by the Josephson current density Jc- The determination 
of the magnetic fleld distribution in a system made by 
two superconducting film forming a planar JTJ is a very 
complex task, since it also involves the Josephson equa- 
tions. The analysis becomes even more difficult when the 
JTJ is biased (the distribution of the bias current peaks 
near the film edges) and the junction dimensions exceed 
the Josephson penetration depth Aj = h/2efj,odeJc- 
However, following Ref. [l4|, we can argue that in the 
independent films regime, the intensity of the in-plane 
magnetic field felt by a planar JTJ placed at the bor- 
ders of two superconducting fllms can exceed by several 
times the value of the applied transverse fleld. In fact, 
the screening currents flowing within a depth A^ ;, in the 
top surface of the bottom film and those flowing within 
a depth A^^t in the bottom surface of the top fllm have 
opposite directions. Consequently the associated mag- 
netic flelds add in the barrier plane and a more efficient 
modulation of the JTJ critical current is achieved. 
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EXPERIMENTS 



In order to prove the advantages to use a transverse 
field rather than an in-plane one, we have measured the 
transverse magnetic diffraction patterns Ic{H±) of pla- 
nar high quality JTJs having different geometries, sizes 
and critical current densities. The samples were placed 
on the axis of a superconducting coil surrounded by a 
Pb shield and a cryoperm can in order to attenuate the 
earth magnetic field. The magnetic field produced by 
the coil was calculated through COMSOL MultiPhysics 
numerical simulations. 
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Overlap type junctions 

Fig. 2 compares the diffraction patterns measured in 
a transverse field of two overlap-type junctions having 
the same geometry and dimensions, but different critical 
current density Jc- The junctions have been fabricated 
with the trilayer technique in which the junction is real- 
ized in the window opened in an Si02 insulating layer. 
The thicknesses of the base, top and wiring layer are 200, 
100 and 400 nrn, respectively. Details of the fabrication 
process can be found in Ref.[16i]. The junction length is 
L — bOOfim, while the width is equal to 4/xm. The base 
and top electrode widths are 540 and 506/xm, respec- 
tively. The so called 'idle region', i.e. the overlapping of 
the wiring layer onto the base electrode is about 3 fim all 
around the barrier. In Fig. 2 the junction critical currents 
have been normalized to their maximum values in order 
to make the comparison easier. The closed circles refer to 
a Nb/Alox/Nb — Nb tunnel junction having a critical cur- 
rent density Jc — QQA/cm?, while the open circles refer 
to a Nb/AlN/NbN - Nb sample with = AOOA/cm'^. 
As expected, considering that the field lines associated to 
the screening current are perpendicular to the electrodes 
edges, the shape of these Ic vs. H± curves looks very 
alike to that expected for long one-dimensional overlap- 
type junctions when a uniform external field is ap- 
plied in the barrier plane in the direction perpendicular 
to the junction length. In fact, according to the analysis 
of Refs. 17, 3|, for small field values (Meissner regime) 



the critical current Ic decreases linearly with the applied 
field IciH\i) = /c(0)(l - \Hi\\/h!), where ij| is the crit- 
ical field at which Ic would vanish if flux quanta did not 
start to enter the junction barrier. The skewness seen 
in the experimental Ic — IIi_ curves (being larger for the 
sample having the larger Jc) is due to the self-field pro- 
duced by the bias current / flowing in a close-by super- 
conducting strip in the chip circuitry. The skewness does 
not prevent us from measuring the junction critical crit- 
ical fields iJ^ determined by linearly extrapolating the 
branches starting at maximum critical current to / = 
(dotted lines in Fig.l); in fact, when / = 0, the bias cur- 



FIG. 2: Magnetic difFraction patterns measured in a trans- 
verse field H± of two overlap-type junctions L x W = 500 x 
i/im^ . The closed circles refer to a Nb/Alox/Nb JTJ hav- 
ing Jc — 60A/cm^, de = ISOnm and \j ~ 50/im. The 
open circles refer to a Nb/AlN/NbN with Jc = 400 A/cm^, 
de — 240nm and Xj ~ 20/im. The dotted lines show the 
procedure used to determine the junction critical field . 



rent self-field effects vanish. The critical field He for a 
long one-dimensional junction in presence of a in-plane 
external field applied perpendicular to the long junction 
dimension L can be expressed as the sum of two terms: 
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The first term in Eq[2] dominates in low Jc sam- 
ples for which L « irXj and corresponds to the criti- 
cal field of a point-like junction which exhibits the well 
known Fraunhofer diffraction pattern. The second term 
H^^ (X \/ Jc/de becomes dominant in the high Jc regime 
when L » ttAj. It was first introduced by Ferrell and 
Prange|l3| in order to describe the self- field limiting ef- 
fects in long inline-type junctions. 

Table. 1 reports the values of Jc and de used to predict 
the critical field h\ from to Eql^l The last two columns 

Sample Jc de L/\j hI = + H^^ 

A/crn^ nm A/m A/m 



Nb/Alox/Nb 
Nb/AlN/NbN 



60 180 
410 240 



10 3.2 
25 13 



36 + 58 = 94 
36 + 131 = 167 



TABLE I: Relevant parameters for two overlap-type Joseph- 
son tunnel junctions whose transverse magnetic diffraction 
patterns are shown in Fig. 2. 
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allow the comparison between 7?^ and He ■ We observe 



that for both samples He > H^ , and the ratio Hc/H^ 
changes from about 30 to about 13 when we move from 
the low to the high Jc junction. The data in Table. 1 un- 
ambiguously indicate that a transverse field can be more 
effective than a in-plane one and, remembering that the 
two samples have the same geometrical details, the effect 
of a transverse field weakens as the junction transparency 
increases. 



R 

mm 



0.08 
0.25 
0.32 
0.50 



R/Xj 



1.6 

- 5 
6.4 
' 10 



A/m 



2.7 
1.5 
1.3 
1.2 



A/m 



18 + 48 = 66 
6 -f 150 = 156 
4-1- 192 = 196 
3 300 = 303 



TABLE II: Relevant parameters of the Nb/Alox/Nb annular 
Josephson tunnel junctions used in the experiments. 



Annular junctions 

Further, we have measured the static properties of sev- 
eral Nb/Alox/Nb — Nb annular JTJs having the same 
critical current density Jc = 60A/cm^ (Aj — 50^m), the 
same width AR = 4^to, but different radius R ranging 
from 80 to SOO/Ltm. The fabrication details are the same 
as those for the overlap-type JTJs discussed previously. 
The diffraction patterns in transverse magnetic field will 
be reported elsewhere j20j; here we focus our interest on 
the values of the critical fields. 

The analogous of Eql2] for a one-dimensional annular 
junction having radius R in an in-plane field is 21|: 



m = H 



H. 



MM 



2.404 



$0 



2lTlJLQdeR 



RJc- (3) 



Again we have a contribution -ff^f , independent on the 
Josephson current density Jc, typical of small and in- 
termediate radius annular JTJs immersed in a uniform 
in-plane magnetic field which results in a periodic radial 
field Hr{9) oc cos6' felt by the junction [23. [23j . In such 
case, the Jc vs. TJy curve follows a Bessel, rather than a 
Fraunhofer, behavior (2.404 is the argument correspond- 
ing to the first minimum of the zero-order Bessel func- 
tion). The second term H^^^'^ = RJc in Eg_pwas numer- 
ically found by Martucciello and Monaco[2l5; considering 
that H^'^' /H^ = ( i^/A.7) V2.404^ it becomes dominant 
when R» Rm^ V2AMXj^^. For given de and Jc, 
Eq|3]has a minimum when R = Rm and linearly increase 
with R when R >> Rm- The last two columns of the 
Table. 2 report, respectively, the transverse critical field 
H^ measured for four annular junctions having different 

radii and the expected parallel critical field ijj according 
to EqOwith Jc = 60A/cto^ and de = ISOnm. 

For all samples we observe once again that ijj > H^ , 
and that the ratio hI / H^ changes from about 25 to 
about 250 when we increase the ring diameter, i.e. the 
top and bottom film widths, confirming that the effect 
of a transverse field strengthens as the electrode widths 
increases. In particular, for the three largest rings, hav- 
ing the so called Lyngby geometry [2J|, i.e. the base and 
top electrode widths match the ring diameter, this ratio 
is proportional to R. [The smallest ring {R/Xj ~ 1.6) 



has the base electrode width equal to 540/im that is con- 
siderably larger than the ring diameter 160/zm.] 



CONCLUDING REMARKS 

We have measured the transverse magnetic diffraction 
pattern of planar JTJs having different geometries, sizes 
and barrier transparencies. Our measurements clearly 
indicate that a magnetic field is more effective to mod- 
ulate the junction critical current Ic when applied per- 
pendicularly (rather than parallel) to the junction plane 
provided the JTJ is fabricated close to the borders of 
superconducting films with a large aspect ratio. This is 
due to screening (or Meissner) currents induced by the 
transverse field that circulate mainly on the film surface 
borders which in turn behave as intrinsic control lines. 
We suggest that a transverse magnetic field can be use- 
fully exploited in those applications where the Josephson 
critical current and the Fiske resonances need to be sup- 
pressed. 
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